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ABSTRACT 

A new mathematical and computational technique for calculating quantum vacuum expectation values of energy 
and momentum densities associated with electromagnetic fields in bounded domains containing inhomogeneous 
^— ^ I media is discussed. This technique is illustrated by calculating the mode contributions to the difference in the 

vacuum force expectation between opposite ends of an inhomogeneous dielectric non-dispersive medium confined 
^1 to a perfectly conducting rigid box. 
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r:^' 1. INTRODUCTION 

Mh| Electromagnetic quantum fluctuation phenomena lie at the heart of many processes where the interface be- 

tween classical and quantum physics plays a prominent role. Among these one may cite areas of quantum 
5 I optics, micro-cavity physics, micro-fluidics, photonic structures, early Universe cosmo-genesis, dark energy and 

^ ' cold-atom technology. In these systems one is often confronted with phenomena that interrelate classical con- 

^^1 tinuum mechanics, classical electromagnetism, cavity quantum electrodynamics and fundamental issues relating 

to fluctuation-dissipation mechanism9^'2l_ jj^ particular, dynamic (material) fluctuations induced by quantum 
fluctuations of the electromagnetic field have experimental consequences and offer an exciting opportunity to 
confront the limitations of basic theory with observable data. In technology such fluctuations may manifest them- 
selves as quantum induced stresses. For example Casimir stresses cannot be ignored as nano-structures develop 
CN i ever smaller miniaturizations. In micro-fluidics, physical processes can be confined to (deformable) micro-cavities 

that are guided by electromagnetic fields. Such micro-laboratories offer new possibilities to explore cavity QED 
experimentally as well as enhancing the control features of micro-fluidic design. Indeed it has even been suggested 
C^ I that chemical processes in such an environment may shed light on the mechanism that evolves inert matter into 

living cells. 



The role of quantum fluctuations in determining the behavior of fabricated micro-structures is becoming 
increasingly important in a wide area of science and technology. Such fluctuations are also at the heart of 
, ■ many fundamental problems in physics ranging from the stability of fundamental constituents of matter to the 

r> I lack of progress in unifying quantum field theory with gravitation. Many of the problems arise due to a lack 

jrt ' of knowledge of basic interactions between fields and matter at some scale and the need to regularize current 

theories in order to make experimental predictions. For renormalizable theories such as QED in the vacuum 
these are remarkably accurate. However macroscopic predictions directly from QED that are affected by the 
presence of bulk macroscopic matter and material interfaces can be made with far less confidence since they 
depend critically on both geometric and constitutive modeling at macroscopic scales^ ■'*. In particular the role of 
quantum states of the electromagnetic field on the detailed behavior of isolated closed material micro-domains 
of polarisable matter remains an unsolved problerrP. 

The quantization of the electromagnetic field in the presence of media satisfying linear electromagnetic consti- 
tutive relations relies on a knowledge of a regular basis of eigen-solutions to a modified Helmholtz type equation 
that determines the electromagnetic fields. If the electromagnetic properties of the media are discontinuous 
across surfaces in space such solutions must satisfy jump conditions across them dictated by Maxwell's equa- 
tions for the fields. The nature of the eigen-spectrum of the vacuum Helmholtz operator on time-parameterized 
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differential forms on space is determined by the global topology of spatial domains with boundaries. Mathe- 
matical procedures exist for analyzing such problems using the Hodge- Weyl-Friedrichs decomposition of forms 
on manifolds with boundary. In principle they can be used to construct a Hilbert space with a real basis of 
transverse (divergence-less) forms (i.e. in the kernel of the co-derivative 6) satisfying Dirichlet and Neumann 
boundary conditions. The split of this space into mutually orthogonal subspaces is responsible for the classifi- 
cation of electromagnetic fields into TE, TM and TEM modes in certain domains. For example, if the vacuum 
of 3-dimensional Euclidean space is partitioned into interior and exterior regions by a closed perfectly conduct- 
ing boundary surface one may establish an ortho- normal basis of transverse 1-forms $^(r), ^J^{r) satisfying 
appropriate boundary conditions in each region. For a hollow closed cavity these basis modes can be defined in 
terms of the eigen-1-forms $Ar(r) and ^Ar(r) of the Hodge-de Rham operator (Laplacian) A on forms in space 
satisfying different boundary conditions: 

A$Ar=/^^$jV, A*Ar = A^*7V, 

where S = —^d^ is the Euclidean exterior co-derivative on 1-forms on a simply-connected domain U G H^, d 
denotes the exterior derivative, —A = d6 + S d and TV labels a triplet of infinitely denumerable discrete numbers 
labelling the real non-zero eigen-values and associated eigen-forms. In a simply-connected domain such 1-forms 
can be employed to represent the Coulomb gauge Maxwell vector-potential 1-form A = A^'^ + A^" where: 

JV JV 

and # denotes the Euclidean Hodge map on forms in space. The eigen-values /xat, \n determine the normal- mode 
frequencies of electromagnetic fields in the cavity. A similar analysis can be performed in the exterior region 
(which may be non-simply connected and involve TEM modes). The computation of the eigenvalues Aat, /i^r 
is the key precursor to all quantum computations since they characterize the extrinsic domain geometry and 
determine the spectral content of the infinite number of quantum oscillators that represent the electromagnetic 
field. However it is often difficult to determine analytically such bases for generic domains. Furthermore if the 
partition involves electrically neutral domains containing linear media that may be inhomogeneous, anisotropic, 
magneto-electric or dispersive this program involves a modified Helmholtz classical boundary value problerrF"^. 

In this paper the modification necessitated by the presence of an inhomogeneous but non-dispersive and 
non-conducting medium contained in a closed rectangular 3-dimensional perfectly conducting cavity is explored. 
This is a precursor to a regularization scheme needed to extract finite quantum expectation values for stresses 
in the medium induced by quantum electromagnetic fluctuations. 

2. INHOMOGENEOUS DIELECTRIC MEDIA 

Consider a smooth open region of space containing a stationary non-dispersive medium characterized by an 
inhomogeneous permittivity scalar e(r) and constant permeability scalar /i = ^o- Denoting time derivatives with 
an over-dot the classical source-free Maxwell system to be solved is: 

de = -B, -db = eE (1) 

M 

dB = 0, d{eE)^0 (2) 

de^O, d/^ = 0, e = 0, /i==0 (3) 

where e, b denote time dependent electric and magnetic 1-forms respectively and E = #e, B = #6. If the 
(time-dependent) spatial 1-form A and spatial 0-form (j) belong to a class of gauge-equivalent potentials defining 
the electric and magnetic fields by 

e = -A-d(j), b = #dA (4) 



then the above system reduces to 



In a particular gauge with 



S{eA) + S{£d(f>) = 
SdA + efi{d(j> + A) = 



<5 (eA) = 
the equation for the scalar potential decouples: 

ddA + e/iA = —efid(j) 
S{ed(j)) = 0. 

Furthermore, we may set d = for systems without free charge^. Hence in this gauge the local Maxwell system 
above is solved in terms of spatial 1-forms satisfying S{eA) — and the equation: 

SdA + efiA = 0. (5) 

This gives rise to a modified Helmholtz equation for time harmonic fields. Across any non-conducting interface 
where the dielectric scalar is discontinuous one has two conditions on e and b restricted to the interface. At each 
point on the interface the jump in the normal component of b and the tangential component of e must vanish. 
Furthermore if there are no real charges or electric currents on the interface the jump in the normal component 
of d where d — ee and the tangential component of h where h — fi~^b must also vanish on the interface. If the 
interface is perfectly conducting one assumes that all fields vanish on the side of the interface that is free of the 
material medium when calculating the jump. It is worth noting that in a bounded inhomogeneous medium where 
e is a continuous function of position one cannot exploit translational invariance in space and spatial Fourier 
transformgi^ to simply express normal modes in terms of eigen-forms of R^ spatial translation operators. 

If a bounded domain U C R^ contains a dielectric with a piecewise inhomogeneous permittivity one writes 
the general real 1-form solution to ([5]) on J7 as 

A(i,r) = ^^^(t,r) (6) 

N 

where N denotes a triple of discrete labels. Suppose the dielectric is composed of M sub-domains where the 
permittivity is em(r) in the ?7i-th sub-domain. Thus 

M 

e{r) = J2 em{r)y„,{r) 

m—l 

where 3^m(T") is unity in the subdomain [/,„ C U and zero elsewhere. For stationary electrically neutral dielectrics 
in domains U bounded by conducting surfaces the electromagnetic jump conditions at interfaces above yield a 
homogeneous system of equations that determine a collection of eigen-modes (up to normalization) and the 
associated eigen- frequencies w^r. The number of distinct eigen-spaces and the degeneracies of the associated 
eigen-frequencies depends on the rank and symmetry of the homogeneous system which in turn reflects how the 
boundary geometry and boundary conditions affect the nature of the global topology of the domain. For a rank 
S system the eigen-modes may be written: 

S M 

s—l m—l 

where the {A}^ ^'^''"(r)} = {^jv (»")}* constitute a basis of solutions to ^ subject to S{eAj^ )^s,m.| _ q ^^^ 

the above jump conditions in the domain U. 



3. ELECTROMAGNETIC ENERGY AND STRESSES 

Classical forces (and torques) transmitted by electromagnetic fields through the vacuum can be encoded into 
a covariant stress-energy-momentum tensor. The modification of such a tensor for fields in a medium has long 
been a subject of debate and experimental investigation. This debate has continued when the fields become 
operators subject to quantum laws. In this article we adopt a symmetrization of the stress-energy-momentum 
tensor for media at rest advocated by Minkowski. In terms of electromagnetic fields in any sub-domain Um the 
instantaneous classical electromagnetic energy is: 

^™= / lieA#d+bA#h). (8) 



The component of the instantaneous integrated electromagnetic stressii", in a direction defined by a unit 
spacelike Killing vector field K that generates spatial translations, transmitted across the side of any portion 
T,„i of an oriented surface in U„i adjacent to the fields in the following integrand, is the force component: 

TK,m ^\l iiK#h A 6 - e A iK#d - #6 A hiK) - d{K) A #e) 

3.1 Computation of Induced Dielectric Stresses by Electromagnetic Mode Fluctuations 
in a Cavity 

The above generalities will now be illustrated for a system comprised of a simply-connected inhomogeneous 
dielectric medium bounded by a perfectly conducting stationary, inextensible rectangular box with sides of 
length Lx, Ly, Lz. Finding non-trivial exact analytic solutions to ([5]) is non-trivial. However, if e(x,y,z) — 
Co P exp{az/Lz) in Cartesian coordinates with real dimensionless positive constants /3,a, then general solutions 
satisfying the above boundary conditions can be expressed in terms of Bessel and trigonometric functions. Since 
the interior U of the box is simply connected the boundary conditions yield S = 2 and a decomposition into 
orthogonal TE and TM modes with respect to the z-axis is possible. With opposite faces of the box at z = 
and z = Lz respectively the TE mode structure is given in the above gauge by: 



-^iv ™(^) = -^^^N^l'nJfiz)] Ikx sm{kxx) cos{kyy)dy - kyCOs{kxx)sm{kyy)dx\ 



where kx — ^^^f^,ky = '^~i ^ stands for the triple (nx,ny,p) with nx,ny positive integers (including zero) and 
A/"™ denotes a normalization constant. Furthermore 



where 






with Cq — -^— in these expressions and the ojjf are the values of the p-th roots of the TE-vaodc spectrum 
generating equation: 

J.];4vWmY^T4vNiLz)]~J^T4vNiLz)]Y^T4vNio)] = o. 

The TAI mode structure is given by: 

^(^+)-™(r) = '^^ ^ '^TWn{^)\ + "" la ( ^- cos(fc,a;) sm{kyy)dx + ky sin(fc,x) cos{kyy)dy 



with normalization constant M™ , 

where 

^^^^ aco ' ^^n) - ^,{k,^ + l^y) + ^. W- ^™(0)r'™[^™(0)] + i;™[r;™(0)]- 

In this case the uj™ are the values of the p-th roots of the TM-mode spectrum generating equation which may 
be written in the form: 

^.™ [^7^(0)] ?.™ [riT{L^)] - Itm [ry^-(L,)] F,™ [r;^-(0)] = 

where f{ri) = rjf'{r}) + /(ry) for any /(ry). These expressions enable one to calculate the field modes for e, 6, h, d 
using ([4]) and the constitutive relations. 

The quantum description can be constructed by generalizing the methods used in vacuum cavity QED. A 
Fock space of quantised modes is introduced by introducing the annihilation and creation operators a^ and a^, 
satisfying the commutator relations: 

[a%, a%, ] ~ 5nn' 5^'^ 

for s £ {TE,TM}. The Fock space vacuum state Ag is annihilated by all aj^. In the above gauge stationary 
quantum modes in a closed cavity are described by the Hcrmitian operator 

M 
se{TE,TM} m=l 

The quantum field modes for the Hermitian operators e, b, h, d follow from (j4]) but with this mode operator and 
the corresponding operator constitutive relations. Replacing the classical fields for the dielectric filled cavity 
in the classical expression for E by such operators yields the quantum hamiltonian £ for quantum fields in the 
cavity. Its expectation value i^A^ [E] in the Fock space vacuum state requires renormalizatiorfi^En] This is effected 
by subtracting from an infinite mode sum an expectation value of the energy of a system with a homogeneous 
medium. To effect this subtraction both mode summations generally require a regularization scheme for their 
definition. Thus for the system with conducting boundaries and a discrete spectrum one defines: 

se{TE,TM} N 

for some suitable smooth function satisfying ■0(0,0;^) = 1 that renders the summations meaningful. Each cavity 
mode labeled by N, s, with eigen- frequency uj%, contributes the factor ^ cuf^ to the vacuum expectation of the 
(regularized) energy. The well-defined condition 

fixes the normalizations Mf^ of the mode amplitudes Aj^ and their conjugates: 



a^^c^ivl^^Y L^Ly IJfnj^ 



where, with fljf — 






IJT = e-{^'l^[ni^e^f-{<^'l^WN]f 



and 



(A/"™)' = 



a^co[{v™Y - 1] L^Ly I™ VI 



where, with f]L*'^ = 2c^^ l.v^ 



(i - (z.^)' + (f^D'e") (<i>?r [f^^e^])' - (i - (^D' + i^N'f) (*r[^? 



The mode contributions to the vacuum expectation values of the induced electromagnetic stress field in the 
dielectric can now be calculated from the Hcrmitian operator-valued stress 2-form: 

^K^\ (iKifi Ab - e AiK#d - #b Ah{K) - d{K) A i^e) . (9) 

The quantum expectation value of the regularized force component in the Fock vacuum state Aq acting perpen- 
dicular to a surface Sz(, of the box a,t z = zq is 



Vzo/reg — -C^Ao 



a a 



rcg 



In a box containing a homogeneous permittivity, the expectation values of the force at the end faces of the box 
aX z — Q,z = Lz are equal. This is not the case for an inhomogeneous dielectric in general. Indeed one finds, 
after some calculation, the difference between the force expectations 

(AJ")rcg = (J'o>rcg - (-T^Ljrcg = ^ X! ^^Ni'i'^^^N) 

^ se{TE,TM} N 

This is a surprisingly simple result given the complexity of the mode structures involved. To effect a renormal- 
ization of this result requires a computation which will not be reported here. 

4. CONCLUSION 

It is expected that a non-zero (AJ^) will survive renormalization when the regulator is removed and this indicates 
that any confined inhomogeneous material dielectric must sustain stresses induced by electromagnetic quantum 
fluctuations if the confining domain is rigid. If the medium remains static such stresses induce mechanical 
(elastic) stresses in the dielectric to maintain equilibrium. Unlike similarly induced classical stresses by the 
classical gravitational field in the laboratory (that vary with the orientation of the dielectric) the quantum 
induced electromagnetic stresses are permanent. In principle they could be detected experimentally by noting 
the variation of the induced stress field within the dielectric with variations of the permittivity inhomogeneities. 
Such variations might be detected using photo-elastic effects on the polarization of light passing through the 
medium. 
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